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ABSTRACT 

This paper is devoted to the analysis of the angular resolution limit (ARL), an important performance measure in the directions-of-arrival 
estimation theory. The main fruit of our endeavor takes the form of an explicit, analytical expression of this resolution limit, w.r.t. the 
angular parameters of interest between two closely spaced point sources in the far-field region. As by-products, closed-form expressions of 
the Cramer-Rao bound have been derived. Finally, with the aid of numerical tools, we confirm the validity of our derivation and provide a 
detailed discussion on several enlightening properties of the ARL revealed by our expression, with an emphasis on the impact of the signal 
correlation. 

Index Terms — Cramer-Rao bound, angular resolution limit. Smith criterion, directions-of-arrival estimation 



1. INTRODUCTION 

As an important topic within the area of signal processing, far-field source localization in sensor array has found wide-ranging application in, 
among others, radar, radio astronomy and wireless communications 1 1]. One common measure to evaluate the performance of this estimation 
problem is the resolvability of closely spaced signals, in terms of their parameters of interest. In this paper we investigate the minimum 
angular separation required under which two far-field point sources can still be correctly resolved. 

To approach this problem, it is necessary to revive the concept of the resolution limit (RL), which will serve as the theoretical cornerstone 
of this paper. The RL is commonly defined as the minimum distance w.r.t. the parameter of interest (e.g., the directions-of-arrival (DOA) 
or the electrical angles, etc.), that allows distinguishing between two closely spaced sources ||2-H]- Till now there exist three approaches to 
describe the RL. The first rests on the analysis of the mean null spectrum (5|], the second on the detection theory yllalZl' and the third on the 
estimation theory, capitalizing on the Cramer-Rao bound (CRB) bBB- A widely accepted criterion based on the third approach, proposed 
by Smith 1 2], states that two source signals are revolvable if the distance between the sources (w.r.t. the parameter of interest) is greater than 
the standard deviation of the distance estimation. In this paper we consider the RL in the Smith's sense, due to the following advantages 
over competing approaches: The Smith criterion /) takes the coupling between the parameters into account and thus is preferable to other 
criteria of the same category, e.g., the one proposed in llsl llOirilll : ii) enjoys generality unlike, e.g., the mean null spectrum approach which 
is designed for a specific high-resolution algorithm; Hi) is closely related to the detection theory approach, as recently revealed in fl. 

This paper investigates the analysis of the RL for two closely spaced correlated deterministic sources. The RL has recently received an 
increasing interest especially after the publication |2] which received the IEEE best paper award. Prior works on the RL consider, on the 
one hand some specific criteria as the RL based on hypothesis tests |3, 12-14], or are tailored to specific estimation procedures as, e.g., the 
MUSIC algorithm in L15.1 . On the other hand prior works based on the Smith criterion either contain non-closed form expressions that require 
numerical evaluation, as in 171 [iql or are based on specific ideal assumptions (e.g., one known DOA I2l. [l7ri . uncorrelated sources llSll , ULA 
case ti|2> J7], non-time-varying sources |P|, etc.) In our work, we propose to derive an analytical expression for the angular resolution limit 
(ARLJj, denoted by 5, between two closely spaced, time-varying (both in amplitude and phase) far-field point sources impinging on non- 
uniform linear array, which, to the best of our knowledge, is till now absent in the cuiTent literature. As by-products, closed-form expressions 
of CRB w.r.t. the relevant parameters are provided to facilitate the derivation of the ARL. Furthermore, our expression, by virtue of its concise 
form highlights the respective effects of various system parameters on the ARL 5. The analysis of expressions for different cases of correlated 
and uncorrelated sources, reveals a number of enlightening properties pertinent to the ARL's behavior, while being also computationally 
efficient, avoiding the difficulties associated with the numerical solution of non-linear equations. 

The following notation will be used throughout this paper: (■)^, (■)"^ denote the conjugate transpose and the transpose of a matrix, 
respectively. tr{-} and vec{-} denote the trace and the vectorization of a matrix, respectively. 5R{-} and £s{-} denote the real and imaginary 
part respectively. 1 1 ■ 1 1 denotes the norm of a vector, ® denotes the Kronecker product, whereas E{-} denotes expectation. 



The so-called ARL characterizes the RL when we consider the angular parameters as the unknown parameters of interest. 



2. MODEL SETUP 

Consider a linear, possibly non-uniform, array comprising M sensors that receives two narrowband time-varying far-field sources S\{t) and 
S2(i), the directions-of-arrival of which are 9\ and 62, respectively. Then the received signal at the m-\h sensor can be expressed as |T|: 

2 
a;™(i) = ^s,(i)e^'=''" ='""**' +n„(t), i=l,...,iV 

and m = 1, . . . , M, 

where the sources are modeled b)OSi(t) — ai{t)e-'^^^''^^^^^ '\ i = 1, 2 in which ai{t) denotes the time- varying non-zero real amplitude, /o 
denotes the carrier frequency, iiiit) denotes the time- varying phase, dm denotes the spacing between the first sensor (which is chosen as the 
so-called reference sensor, i.e., di = 0) and the ra-th sensor, fc = ^ is the wave number (with A denoting the wave length), rim (t) denotes 
the additive noise at the m-th sensor, and A^ is the number of snapshots. 

Fore mathematical convenience, we define Ui — k sm{6i), i = 1, 2 as our parameters of interest. Changing ([I) into the vector form, one 
obtains: 

x{t)^As{t) + n{t), (2) 

where a; (f) — [xi{t), . . . ,XM{t)] , s{t) = [si{t) , S2{t)] , n{t) — [ni{t), . . . ,nM{t)] , and A — [a{i^i),a{i'2)]- The steering vectors 
are defined as a{i'i) = [e-''^''^^ , ■ ■ ■ , e-'"''^*^]"^, i = 1, 2. Furthermore, define the correlation factor p between the two signals as L20,1 

P = "n — in — ii (3) 

l|si|| I|S2|1 

where Si = [sj(l), . . . , Si(A'^)] , i = 1, 2 are signal vectors. 

The following assumptions are made in the remaining of this chapter: 

Al The sensor noise follows a complex circular white Gaussian distributed, both spatially and temporally, with zero-mean and unknown 
noise variance a^. 

A2 The source signals are assumed to be deterministic; and the separation of the sources is small. 

A3 The unknown parameter vector is ^ = [i/i , 1/2 , cr J . Thus, for given ^, the joint probability density function of the observation 

X = [a;^(l),...,a;^(Af)]^ can be written as p(x | £.) = ^mn^h^ exp (-(x - f£)"R'^(x - m)). where .R = o'^ImnwiA 

/.= [(As(i))^,...,(As(iv)ri^. 



3. DERIVATION OF 5 

The derivation of the ARL 5 can be divided into three steps. The first step involves the derivation of the CRBs w.r.t. the relevant parameters. 
The second builds on this result and simplifies the implicit function based on the Smith criterion, the root of which yields 5. The last step is 
to solve the function corresponding to different values of p, leading to the final expression for the ARL. 

3.1. CRB Derivation 

The CRB of the unknown parameters {ui and V2) is obtained as the analytical inverse of the Fisher information matrix (FIM) for ^ (denoted 
by I). Under Gaussian noise, the elements of X can be calculated using the following expression |21]: 

where [^]^ denotes the i-th element of the parameter vector ^. Thus for our model, I takes the following block-diagonal form: 



where 
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(6) 



and 



, which a = Y.m=i dL, SNR, = e?/c^^ i = 1, 2, where e^ = sJy^I^oW)/^ ^ i = 1, 2; 1 

M M 

V = sf S2 ^ d^e"-'''"''^^^'^^' = s"s2 ^ d^e"-'''™^, (7) 



-Note that this is a commonly used signal model in communication systems (cf. iTsLfT^I ). 



where A = i^i — 1^2 denotes the spacing between i/i and 1/2- We assume in the following that i/i > U2, hence A > 0. 
By inverting the 2x2 matrix I we obtain the following expressions for the entries of the CRB matrix: 

2Na. 



2Na 
UKli(!/2) = {'2. I = 
L J2,2 

and 



CRB(/.i) ^ [t-'] ^ ^ = ^SNR2, (8) 

CRB(/.2) ^ \t-'] = ^SNRi, (9) 

L J 2,2 W 



CRB(j.i,;.2)^ [X ']^^ = -^5^5R{77}, (10) 

where * = 4q27V^SNRi ■ SNR2 - (4/^*) ■ ^'^{v} is the determinant of I. 

3.2. Equating the ARL 

According to the Smith criterion, the ARL, S, is given as the angula r spacing. A, which is equal to the standard deviation of the estimate of 
A. The latter, under mild conditions ||22|| . can be approximated as \/CRB(A), suggesting that S can be obtained as the (positive) solution of 
the equation: 

S^ = CKB{5). (11) 

where CRB{S) = CRB(i/i) + CRB{u2) - 2CRB(!/i, 2^2) fl- 
Substituting iSt-ilOb into dllb . the latter is transformed into: 

5^ = CRB(i^i) +CRB(i/2) - 2CRB(i/i,zy2) 

2 2 (12) 

= -{N ■ SNR2a + N ■ SNRia + — 3ff{r?}). 

Substituting S for A in identity (O we observe that ( I12t is a highly non-linear equation in 5. Hence, in order to obtain the solution of | |121 
w.r.t. 5, and taking into account that S is small, we resort to the first-order Taylor expansion of rj around (5 = (Qto obtain: 



ri ^ sf S2 ^ dm(l - jdmS) 

m — 1 
/ M M \ 

= sf S2 i^dl,-j5^di^\ 



(13) 



(14) 



= sfs2(Q-i5/3), 
where /3 = J2^=i d^i- Combining l[T3) with (O, it follows that: 

5RW~l|si||||s2||K{p(Q-j5/?)} 

= A''eie2(pQ + p/35), 

where p and p are defined as the real and imaginary part of p, respectively, i.e., p = 5R{p} and p = T{p}. 

Now we merge ( 114b into ( 112b and, after some mathematical manipulations, obtain the following quartic function of 5: 

D^5^ + 2CD5^ + {C' - AB)5' + D5 + ^-^ -^ C = 0, (15) 

where A, B,C and D are defined as: 

A = N ■ SNRia, (16) 

B = N ■ SNR2Q, (17) 



C = A^VSNRi ■ SNR2pa, (18) 

and 

D = TVVSNRi ■ SNR2P/3. (19) 

Thus our task of finding the expression of S has been brought down to finding the root of ( 115b . 



^In asymptotic cases <5 becomes very small and our approximation made here is tight, as will be proved by our simulation (cf. Fig.[T]. This can be explained 
by the fact that the Maximum Likelihood estimator, and generally all high resolution estimators, have asymptotically an infinite resolution capability leading 
to <5 -i> EHH. 



3.3. Expression of 5 for different correlation factors 

The solution of dlSl l. depending on different values of p, falls into the following three cases: 

Case 1. Non-zero imaginary part of the correlation coefficient p(p 7^ Q): in this case M5\ remains a quartic function in 5. We know from the 
parameter transformation property of the CRB (cf [Q, p.37) that CRB(5) = CRB(-5), Thus, if (5 is a root of ([TT) (hence of ([15)), 
then ~5 will also be a root thereof, which allows us to remove all the terms of odd degrees in ( 115) . leading to a quadratic equation of 
5^ . Hence 

— + C = 0. (20) 



D^J^ + {C^ 



AB)5^ + — 



The root of (Hoj ifl 



AB 



C^ ~ ^{C^ - AB)2 



4D2(il±s +(7) 



2D2 



^(1- 



.,2 



(21) 



where 7 = (1 - p^)a^, k = 2/9^/3^ and 



+ 



+ 



2p 



AT VSNRi SNR2 VSNRi ■ SNR2 



(22) 



The existence of 5^ in ( I21t is assured since under realistic conditions {aKcfi/'y^) <C 1. Thus the ARL is given by: 



\ 



anq, 

_,2 



(23) 



Case 2. Not fully correlated signals with zero imaginary part of the correlation coefficient p (p = and p 7^ ±1): in this case D = 0, 
(C^ - AB) ^ 0, and ([BJ degenerates into {C^ - yl_B)(5^ + ^^ + C = 0. Taking its positive root we have: 



(5 = 



2 



+ C 



AB-C2 



(24) 



The existence of 5 is guaranteed from the fact that in this case both and 7 are greater than zero. It is worth noticing that an important 
special case of Case 2, in which both p and p equal zero, namely, the two signals are uncorrelated, reduces ( 124b to 



1 



1 



2Na VSNRi SNR2 



(25) 



Case 3. Fully correlated signals with zero imaginary part of the correlation coefficient p (jj — and p — H): in this case ( |15t degenerates 
to —5-^ + C — and a solution can not be found|j 

Now, combining the results of all three cases presented above yields our final expression of the ARL that can be written as: 



5={ 




f or p / 



f or p = and p 7^ ±1 



(no closed-form expression available), f or p = and p = ±1 



(26) 



Note that, for the uniform linear array (ULA) configuration the parameters in 

71 ,r2 / 71 //" -1 \ 2 

-^ — —d, where d denotes the inter-sensor spacing. 



can be derived as a 



M(M-1)(2M-1) 
6 



d and /3 



"^The other root of ( I20K which is very large, is in contradiction with the observation made in Footnote 3, thus is regarded as a tiivial solution and rejected. 

'One can expect that for the case in which p = ±1, i.e., the two signals are linearly dependent, the approximation made using a first order Taylor expansion 
is not sufficiently tight w.r.t. the true model. Thus in this case it entails a higher order Taylor expansion and thereby involves solving a sextic equation, the 
detailed analysis of which, unfortunately, is due to the space limitation beyond the scope of this paper 



4. SIMULATIONS AND NUMERICAL ANALYSIS 



The context of our simulation is a ULA of M 
N = 100. Our results are as follows: 



6 sensors with half-wave length inter-element spacing. The snapshot number is given by 



In Fig.[T]we validate our approximate analytical expression of 5 in l |26t for two cases (p ^ 0; p — O&p 7^ 1) by comparing it with the 
true S (obtained by solving Jilt numerically) and show that both results are identical. 



As is revealed by ( I26t . the concrete waveforms of the signals has no effect on S, which only depends on the two signal's respective 
strengths (ei, £2) and the correlation p between them. Furthermore, note that either p or p plays its role separately. Fig.|2]shows that, 
with a fixed p, the ARL 5 slightly increases with the value of |p| . However, this impact is so limited compared to that of the parameter 
p, such that the former is practically negligible, (cf. Fig. |3] and Fig. |4] both of which show that 5 increases notably as p raises, while 
remains nearly unaltered with the change of p.) This fact can be explained by considering, for (aKcp/j'^) <C 1, the first order Taylor 
expansion to l |23l l around [ancf)/')'^) = that is given by: 



!(,_(, _^ 



(27) 



which is the same expression as J24t . independent of p. 

The dependence of 5 on the signal strengths is reflected in the expression of (cf. ||22I(). where we see that if the strength of one 
signal is much greater than the other, e.g., ei 3> £2, then (fi ~ i/{N ■ SNR2), and S becomes restricted by the weaker signal. Thus 
enhancing the strength of only one signal cannot infinitely diminish 5, as is shown by Fig. [5] in which we increase ei from 1 to 1000 
while keeping £2 = 1, and find that S converges to a certain value (determined by £2). 

Fig. |5] also investigates the impact of the sensor array geometry on 5 (cf. Table 1) and reveals that a loss of sensors in the array 
configuration has a considerable impact on 5 only when it causes a diminution of the aperture size of the array, as in the case of Type 
1. If, however, the array aperture remains unchanged, as in the case of Type 2, this impact is considerably mitigated. 



Array Type Geometric Configuration 



Type 1 
Type 2 
Type 3 



o»»o»»oo 



• oooo»»» 



Table 1. Different aiTay geometric configurations. • and o represent the position of sensor and missing sensors, respectively. The inter- 
element spacing is half-wave length. 
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Fig. 1. Numerical and analytical & vs. a for £i = £2 = 1, with p = 0.5 + 0.5j and p = —0.5 respectively. 



5. SUMMARY 



In this paper we studied the angular resolution between two point sources and provide a closed-form expression for the ARL 5, the validity 
of which was confirmed by simulation. We also noticed that 5 is not dependent on the special waveforms of the signals, but only on their 
strengths and the correlation factor between them, and that the imaginary part of p only has a negligible impact on 5, while the impact of the 
real part of p is decisive. Furthermore, we showed that 5 is constrained by the weaker signal, and therefore cannot be infinitely decreased. 
Finally the impact of different array geometries on & is discussed. 
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Fig. 2. 5 vs. \p\ for p = 0.5, £i = £2 = 1 and a^ = 1. 
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Fig. 3. 5 vs. \p\ for £i = 62 ~ 1, o-'^ ~ 1 and various p. 
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Fig. 4. S vs. p for ei = £2 = 1, cr^ = 1 and various p. 
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Fig. 5. S vs. £2 for p — 0.5 + 0.5j, £1 = 1, a =1 and various array configurations. 
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